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FLOCCULENT FLOWS IN THE CHROMOSPHERIC CANOPY OF A SUNSPOT 
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ABSTRACT 

High-quality imaging spectroscopy in the Ha line, obtained with the CRisp Imaging SpectroPolarimeter 
(CRISP) at the Swedish 1-m Solar Telescope (SST) at La Palma and covering a small sunspot and its surround- 
ings, are studied. They exhibit ubiquitous flows both along fibrils making up the chromospheric canopy away 
from the spot and in the superpenumbra. We term these flows "flocculent" to describe their intermittent charac- 
ter, that is morphologically reminiscent of coronal rain. The flocculent flows are investigated further in order to 
determine their dynamic and morphological properties. For the measurement of their characteristic velocities, 
accelerations and sizes, we employ a new versatile analysis tool, the CRisp SPectral EXplorer (CRISPEX), 
which we describe in detail. Absolute velocities on the order of 7. 2-82. 4km s _I are found, with an average 
value of 36.5±5.9km s _1 and slightly higher typical velocities for features moving towards the sunspot than 
away. These velocities are much higher than those determined from the shift of the line core, which shows 
patches around the sunspot with velocity enhancements of up to 10-15 km s _1 (both red- and blueshifted). Ac- 
celerations are determined for a subsample of features, that show clear accelerating or decelerating behavior, 
yielding an average of 270+63 m s~ 2 and 149±63 m s~ 2 for accelerating and decelerating features, respectively. 
Typical flocculent features measure 627+44 km in length and 304+30 km in width. On average 68 features are 
detected per minute, with an average lifetime of 67.7+8.8 s. The dynamics and phenomenology of the floc- 
culent flows suggest they may be driven by a siphon flow, where the flocculence could arise from a density 
perturbation close to one of the footpoints or along the loop structure. 

Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: activity — Sun: sunspots 



1. INTRODUCTION 

The advent of imaging spectroscopy at high-resolution 
solar telescopes using adaptive optics and image post- 
processing allows the study of solar phenomena at unprece- 
dented spatial, temporal and spectral resolution. In this pa- 
per we analyze a high-quality set sampling the chromosphere 
above a small sunspot and its surroundings in the Ha line. 
The chromospheric canopy exhibits ubiquitous fast-moving 
small-scale features constituting flows both towards and away 
from the sunspot in the superpenumbra, as well as in fibrils 
anchored in plage only. We call these features and flows "floc- 
culent'Q to characterize their intrinsic intermittency. 

Similar blo b-like phenomena have recently be en ob - 
served by [Van Noort & Rouppe van der Voortl (120 06), 
ISanchez-Andrade Nuno et all (120081) . IWedemeyer-Bohrnl 
( 120101) . and iLin et al.l ( 120111) . and have been interpreted 
as a manifestation of propagating waves rather than actual 
mass motion. With spatial dimensions of 0.5-1 arcsec or 
less, the features observed are seen to propagate along the 
fibrillar canopy structure observed in Ha and Ca n 8542 A, in 
some cases with additional periodic lateral displacement, at 
velocities close to (and in most cases well in excess of) the 
chromospheric sound speed. 

The flocculent flows above the sunspot and in its super- 
penumbra that are reported on here are reminiscent of coronal 
rain as well as the inverse Evershed effect, and may consti- 
tute mass motions rather than waves. We therefore complete 
this introduction with brief reviews of the coronal rain phe- 
nomenon and the (inverse) Evershed effect. 
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1 They are not the "flocculi" reported in th e older solar literature, in refer- 
ence to unresolved Ca n H & K network (e.g. Bray & Loughhead 1974). 



Coronal r ain — In coronal loops both (transient) siphon 
flows (e.g iKieldseth-Moe & Brekkejl99llDovle et alJ2006t 
iTian et all 120091) and "blobs" of cooled down plasma (so- 
called cor onal rain) precipitating along the legs of coronal 
loops ( e.g..lSclu-iiverf200lt|Muller et alj2005UDe Groof et all 
120051: lO'SheaetalJ 120071) have been observed. Typi- 
cal velocities for coronal rain range bet ween a few tens 
to a few hundred ki l ometers per second dSchriiveri 120011 : 
iDe Groof et ail 120051: IP' Shea et al.l 12007b. wh ich was cor- 
roborated throu gh simulations by Mtiller et al.] (20051). Re- 
cent studies by lAntolin & Rouppe van der Voort 120121) and 
lAntolin et alj ( 120121) have further expanded the statistics for 
both off-limb and on-disk cases, respectively, confirming ear- 
lier dynamics results. Catastrophic cooling is generally ac- 
cepted as explanation of the coronal rain ph enomenon. In 
this process, as advanced by ISchriiveri d2001l) . thermal non- 
equilibrium leads to condensation of plasma at the loop apex 
followed by a runaway cooling process, with subsequent 
downflow of one or more blobs along one or both loop legs. 
The siphon flow velocities in coron al loops are similar to the 
velocities found for coronal rain (cf. lKieldsefh-Moe & Brekkel 
119981: iDovle et al.l 120061) . but do not result in blobs falling 
down towards the photosphere. 

(Inverse) Evershed effect — The inverse Evershed effect is 
a chromospheric inflow of plasma along superpenumbral 
fibrils towards the sunspot umbra, while the outflow ob- 
served in the photospheric penumbra is referred to as the 
(normal) Evershed effect . The latter was first discov- 
ered by lEvershed dl909al) and has been extensively stud- 
ied in the years since, leading to relatively well estab- 
lished values for the outflow velocity of the order of a 
few kilometers per second near the umbra-penumbra bound- 
ary increasing up to about 10 km s at the outer penum- 
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bral edge. Though several models have been proposed in 
the past, thermal convective motions are nowadays most 
widely accepted as the driving mechani s m for the Ever- 
shed e ffect (cf. ISpruit & Scharmeri 120061: iHeinemann et all 
2007t iRempel et alJ l2009t iNordlund & Scharmer! 120101: 
Scharmeretal]l2011l) . 

The existence of the inver s e Ever shed effect was al- 
ready suggested by lEvershed (1909b) following his dis- 
covery of the photospheri c outflow and more extensively 
studied by ISt. John! (119131) . who confirmed the earlier find- 
ings of Evershed that in the chromosphere the flow re- 
verses sign with respect to the photosphere. Siphon flows 
are most often proposed a s the mechani s m behind the 
inverse Evershed effect (cf. IHaugenl [l969t iMaltbvl [1971 
Alissandrakis et alj|1988t iDere et al.lll990t iBoerner & Kneer i 
1992[ iKjeldseth-Moe et alJ Il993t iMontesinos & Thomas! 
199% . with only few alte rnatives, such as gravitationally 
driven in- and downflow dBeckersI 119621) or some mecha- 
nism analogous to the moving flux tube model dTeriaca et al. 
120081). A wide range in velocities is found, however. ISt. John 
dl913l) found the inverse Evershed effect to propagate with 
about 3 k m s" 1 on average and si milar velocities were ob- 
tained by iBones & Maltbvl (119781), whil e somewhat higher 
velocities were found b v IHaugenl (fl%9l bere et alJ (fl990h 
and lAlissandrakis et all d 1988b . However, measurements of 
single a b sorbin g elements in the Ha spectra analyzed by 
Haugen (1969) resulted in velocities of up to 50km s _1 
and ev en higher velocit ies w ere obtained in filtergram stud- 
ies by iBeckersI (fl962h and IMaltbvl (fl975h . This differ- 
ence between spectroscopic an d filtergram velocit i es wa s 
already pointed out before by [Alissandrakis etld] dl988l) . 
The inverse Evershed effect has also been observed in the 
transition region, albeit at typically higher Doppler veloc- 
ities on the orde r of a few tens of km s' 1 obtained from 
spectral analyses d Alissandrakis et al.ll 1 988t fPere et al.ll 19901 : 
IKjeldseth-Moe et alJll993l:lTeriaca et alJl2008l) . and it has not 
been as extensively studied as in the chromosphere. 

In summary, the basic picture that arises from previous 
studies of the (inverse) Evershed effect is one of relatively 
slow penumbral outflow in the photosphere, a faster inflow 
in the superpenumbra at chromospheric levels, and an even 
faster inflow in the transition region. 

In this paper we describe, measure, and interpret the 
flocculent flows in our data. The measurements on the data 
were done using a newly developed analysis tool, the CRisp 
SPectral EXplorer (CRISPEX, further details in Appendix©, 
that is primarily aimed at imaging spectroscopy data. The 
observational data are introduced in the following Sect. |2l 
after which the results of a statistical analysis of these flows 
are presented in Sect. [3] Section [4] offers a discussion of the 
results obtained in the light of similar and possibly related 
phenomena such as the inverse Evershed effect and coronal 
rain, after which conclusions are drawn in Sect. [5] 

2. OBSERVATIONS, DATA REDUCTION AND ANALYSIS METHODS 

2.1. Observational Setup 

The Ha data used in this study were ob tained with the 
CRisp Imaging SpectroPolarimeter (CRISP, IScharmer et alJ 
20081). at the Swe dish 1-m Solar Telescope (SST, 
Scharmer et alJ l2003al) at La Palma (Spain). In the setup 
for these observations the light from the telescope is first 
guided through an optical chopper and a wavelength selec- 



tion prefilter before entering the CRISP instrument (a dual 
Fabry-Perot interferometer (FPI) system with a transmission 
FWHM at the Ha wavelength of 6.6 pm and wavelength 
tuning within a spectral line on the order of <50ms). The 
chopper ensures synchronization of the exposures obtained 
and for the observations considered here, the Ha line was 
selected by using a prefilter centered on 6563.8 A with a 
FWHM of 4.9 A. After CRISP, the light beam is split by 
an orthogonally polarizing beam splitter onto two cameras. 
In addition, before entering the FPI but after the prefilter, 
a few percent of the light is redirected to a third camera, 
which serves as a wide-band anchor channel in the image 
post-processing. All three cameras are high-speed low-noise 
Sarnov CAM 1 Ml 00 CCD cameras operating at a frame rate 
of 35 frames per second with an exposure time of 17 ms. 

2.2. Data Acquisition and Reduction 

The Ha sequence in question was acquired on 2008 June 
1 1 between 7:55-8:32 UT. The data consist of a time series of 
profile scans sampling the Ha line at 23 wavelengths, ranging 
between +1.1 A around line center at 0.1 A spacing. At each 
wavelength a burst of 8 exposures was taken, thus resulting 
in an overall cadence of 6.2 s per profile scan. The field-of- 
view covers an area of about 67x67 arcsec 2 and contains a 
small sunspot with a rudimentary penumbra (AR 10998), lo- 
cated at (X, Y) x (-684,-152), corresponding to p-0.67. Fig- 
ure Q] shows samples of the observed region in Ha line center 
{left panel) and the blue wing (Ha-0.7 A, right panel). The 
pixel size is 0.071 arcsec px -1 , well below the SST's Rayleigh 
diffraction limit for Ha, which is 0.17 arcsec. 

With the use of real-time t ip-tilt correction and t he adap- 
tive optics system at the SST (Scharmer et al. 2003b), as well 
as the image post-processing technique Mul ti-Obi ect Multi- 
Fram e Blind Deconvolution (MOMFBD, iVan Noort et alJ 
120051) . the quality of observations can be greatly improved, 
the latter being employed in order to remove the remaining 
high-order seeing effects. To that end, all images (i.e., at 
each wavelength within the profile scan) were divided into 
64x64 px 2 overlapping subfields to be processed as a single 
MOMFBD restoration. In that process, the wide-band ex- 
posures act as anchors that enable the precise alignment of 
the restored narrow-band CRISP images. More information 
on the MOMFBD post-processing of similar data s ets can be 
found in I Van Noort & Rouppe van der Voortl d2008l) . 

After MOMFBD restoration, the data was further corrected 
for the transmission profile of the wide-band prefilter and for 
the effects of diurnal rotation of the image that result from the 
alt-azimuth design of the SST. F i nally, the images are also de- 
stretched following IShine et alJ d 19941) . which removes most 
of the remaining small-scale rubber-sheet seeing effects. 

2.3. Flocculent Flows 

Closer inspection of the image sequence in either wing of 
Ha shows the presence of flows along what, by comparison 
with the line center image, appears to constitute the chromo- 
spheric canopy (including as such the superpenumbral fibrils 
of the small sunspot). The flows are neither uniform nor con- 
tinuous but rather flocculent, i.e., the flows are observed as 
dark condensations or blobs flowing along the chromospheric 
fibril structures. 

Figure [2] (and in particular the movie of that figure) shows 
examples of flocculent features (arrows). Note that the fea- 
tures in the middle left panel are not visible in the right- 
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Fig. 1. — Sample images. Left: full field-of-view in Ho- li ne cen ter. The loop paths defined with CRISPEX along which space-time diagrams have been 
extracte d are overlaid. The dashed path is discussed in detail in !3. 1.1 l and Fig. [4] The numbered white boxes indicate regions of interest that are discussed further 
in Sect. I2.3l (cf. Figs. [2] and [3] respectively). The arrow at the top of the panel indicates the direction of the limb and the line-of-sight (pointing away from the 
observer). Right: same field-of-view in Ha— 0.7 A. 
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Fig. 2. — Cutout from Fig. [T] of region of interest (ROI) 1 at different wavelengths: Hff-l.lA (Av = -50.2km s , leftmost panel), Ha -0.7 A (Av = 
-32.0km s , middle left panel), Ha line center itself (middle right panel) and Ha+0.7A (Av = +32. 0km s , rightmost panel). The arrows indicate some 
examples of inflowing and outflowing features. A movie of this figure is provided electronically. 



most panel and vice versa, i.e., the same features are not 
visible in both wings simultaneously. The apparent coinci- 
dence of the trajectories of these flocculent features with su- 
perpenumbral fibrils as well as the main flow direction are 
reminiscent of the properties of the inverse Evershed flow. 
It should be noted, however, that these flocculent flows oc- 
cur more generally (i.e., also along fibrils that connect dif- 
ferent network locations) and not necessarily only propagat- 
ing towards the sunspot. Although some of these outward 
moving features can be found near the sunspot, they typically 
appear further away from the sunspot, towards a downflow 
footpoint anchored in the network. The features occur typi- 



cally in streams, either throughout the whole observing period 
or, what appears to be more common, in recurrent, shorter 
episodes of intermittent appearance. 

Some fibrils remain dark in the wings of Ha throughout 
most of the time series or show slower dynamics than the floc- 
culent flows. The right-hand panel of Fig. Q]shows a number 
of these. In some cases the flocculent features and these con- 
stantly present dark fibril structures seem superimposed on 
each other, whereas in other cases also fibrils move down to- 
wards the sunspot, albeit at smaller velocities (cf. for instance 
the middle left panel of the online movie). 

Figure [3] shows spectral properties of a typical flocculent 
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feature. The bottom panel shows the feature in an intensity 

Doppler velocity [km s '] 
-20 20 




1200 
















1000 












\ \> 




800 




V \ 


rr 










600 
















400 








200 




/ N 








^ s 















r ^. 








6562.0 6562.5 


6563.0 6563.5 



Wavelength [A] 




Fig. 3. — Spectral properties of a typical flocculent feature occurring in 
ROI 2. Top panel: spectrum- time diagram at the location marked with a 
white cross in the bottom panel, t = is defined as the instant sampled by 
the middle and lower panels. Around ( = —50 s another flocculent feature 
crossed that location. The white dashed line indicates the average line core 
shift over the displayed time range. Middle panel: spectral profile (solid line) 
and difference profile (the spectral profile subtracted from the field-of-view 
averaged profile, dashed line) of the flocculent feature. A similar set of pro- 
files {dash-dotted and dotted lines, respectively) is shown for the neighboring 
location, without a flocculent feature, specified by the cross-diamond in the 
bottom panel. In addition, the the full field-of-view averaged profile (solid 
curve with diamonds) is shown. Bottom panel: intensity image at Ho- -0.7 A, 
showing an example of a flocculent feature. The dashed curve specifies the 
path of the flocculent feature. 

image in the wing of Ha at a shift of -0.7 A (or equivalently, 
Av = -32.0km s~'). During its visibility in the data sequence, 
the feature moves along an arched path from the left to the 
right through the field-of-view, curving down to the lower 



right of the subfield shown (towards the sunspot which lies 
just outside this subfield on the lower right side), with similar 
other features preceding it and features in adjacent fibril struc- 
tures. The middle and upper panels of Fig.[3]illustrate that the 
flocculent features are most clearly observed in the Ha wings. 
The flocculent flows are composed of features with typical 
Doppler shifts between 20-40 km s in either the blue or red 
wing of Ha, depending on their propagation direction with 
respect to the observer. Note that, the features appear to show 
up as a separate component to the local spectrum (the dashed 
line), as illustrated in the middle panel. In this particular case, 
where an inflow towards the sunspot is observed, the separate 
component appears to be superimposed upon an overall shift 
of the line core. Though the extent to which this separate com- 
ponent is visible in the spectrum may vary with location, time 
and propagation direction with respect to the line-of-sight, it 
is a typical property of the flocculent flows. 

2.4. Methods 
2.4. 1 . Velocities and Accelerations 

Using the path-drawing functionality in CRISPEX, floccu- 
lent features were manually traced along their visible trajec- 
tories and space-time diagrams were subsequently extracted 
along these paths. The left-hand panel of Fig. Q]shows the fi- 
nal selection of the determined paths overlain on the Ha line 
center image, while Fig.|4]shows such an extracted space-time 
diagram corresponding to the dashed track in the former fig- 
ure. As the features appear as dark streaks in these diagrams, 
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Fig. 4. — Part of the space-time diagram extracted along dashed track 
in Fig. [T] The space-time diagram was extracted at Ha ~0.7 A (Av = 
-32.0km s ). The dashed curves indicate measurements of the projected 
velo city a nd acceleration components and illustrate the method discussed in 
Seet.lZ3l 

the determination of their projected velocities and accelera- 
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tions is a relatively trivial task, that was further simplified by 
using the CRISPEX auxiliary program TANAT. 

The majority of the dark streaks in the extracted space- 
time diagrams are straight, and hence a single measurement 
suffices to obtain the projected velocity component in those 
cases. However, some features have curved space-time tracks; 
these are fitted with a parabola to three points along them, 
yielding estimates of the corresponding acceleration or decel- 
eration. Dashed lines have been overlain in Fig.|4]to illustrate 
this procedure. A measure for the error in the projected ve- 
locities and accelerations was determined by repeatedly mea- 
suring randomly-selected dark streaks in each space-time di- 
agram resulting in a measurement uncertainty of 4.0 km s _1 
and 63 m s~ 2 for the velocity and acceleration, respectively. 

The Doppler velocity of the separate component was ob- 
tained employing two methods, a determination of the first 
moment with respect to wavelength and a Gaussian fit, both 
to the difference pro file (cf. the middle panel of Fig. |3] >. In 
the first we followed iRouppe van der Voort et al.l (120091) (see 
also Eq. (1) therein), where the first moment method was ap- 
plied to Rapid Blue Excursions (RBEs). The integration is 
only performed over the range where I™ 8 - lx is positive (i.e., 
the feature is in absorption). For the Gaussian fit method, the 
same constraint is enforced and only the spectral points over 
that same range are used in fitting a single Gaussian with the 
height, position and width of the peak as the three parame- 
ters to be fitted. For both methods, an estimate of the error in 
each Doppler velocity determination was obtained by repeat- 
ing the calculations with boundaries progressively shifted in- 
wards by one wavelength position at a time. The error for that 
particular Doppler velocity was then set to the standard devi- 
ation of those measurements, resulting in an average error of 
0.7 km s and 4.3 km s _1 for the first moment and Gaussian 
fit method, respectively. 

While the separate component constitutes in large part the 
line-of-sight signal of a flocculent feature, it is typically su- 
perimposed on an overall shift of the line core, which may be 
indicative of a lower velocity flow at the same location. Both 
a Gaussian and a 4 th order polynomial were therefore fitted 
to the core of the line (typically the middle 11 wavelength 
positions) at each pixel and for each time step. Comparison 
of the resulting fits with the actual local line profile showed 
the 4 th order polynomial to be more robust in fitting the core 
of the Ho- line than the Gaussian. Hence, the results from 
the polynomial fits were used to construct a Doppler map for 
each timestep, where the offset between the local minimum in 
the line profile and that in the average line profile was con- 
verted to a Doppler shift velocity. This is also a habitual way 
in which velocities for the inverse Evershed effect have been 
obtained in spectroscopic studies. 

2.4.2. Morphological Parameters 

The flocculent features' lengths and widths were deter- 
mined as follows. For each flocculent feature, the intensity 
profile along the spatial axis was extracted from the space- 
time diagram at the wavelength closest to its average Doppler 
velocity and subtracted from the field-of-view averaged inten- 
sity at that wavelength. Next, a difference profile was deter- 
mined by subtracting the average background intensity pro- 
file (i.e., the profile averaged over the time intervals in the 
space-time diagram where no flocculent features were mea- 
sured) from the feature's intensity profile. The difference 
profile was subsequently lowered by the smallest of the two 



minima closest to the maximum representing the feature lo- 
cation, as not doing so typically resulted in an overestimation 
of the FWHM. Finally, a Gaussian fit to the resulting profile 
was performed, yielding a FWHM that was converted to a 
length. However, in order to ensure that only flocculent fea- 
tures (rather than the background) were measured, those mea- 
surements for which the feature FWHM was smaller than the 
FWHM resulting from a fit to the background profile (over 
the same spatial points) were rejected. This procedure was 
then repeated for each point along the velocity measurement, 
yielding a set of length measurements for each flocculent fea- 
ture, the average of which was taken to be the length of the 
feature in question. A similar procedure was adopted to ob- 
tain feature widths, but using the intensity image instead, i.e., 
stepping along the traced path of the feature while it prop- 
agates and extracting the intensity profile perpendicular to 
(rather than along) the propagation direction. 

Since the values of the fitted parameters are influenced by 
the number of points included in the Gaussian fit (i.e., the 
more points included, the larger the resulting FWHM typi- 
cally is), both the length and width profile were fitted over 
three different ranges. Lengths were calculated by taking 11, 
2 1 and 3 1 points into account, while the width profile was fit- 
ted with 5, 11 and 17 points, in both cases centered on the 
coordinates where the feature was measured. The FWHM re- 
sulting from the different fits for which the 1-cr value is lowest 
and equal to or smaller than 1 px was then selected as the fea- 
ture size, while the standard deviation between the three re- 
sults was taken to be a measure of the uncertainty in the size 
determination. The average error for the lengths is 44 km, 
while that for the widths is 30 km. 

3. RESULTS 
3.1. Dynamics of the Flocculent Flows 

3.1.1. Projected Velocities 

In total we obtained 259 velocity measurements of which 
188 can be identified with features moving towards the spot, 
the remaining 71 features moving in the opposite direction. 
As pointed out before, the latter features are observed mostly 
towards the footpoints in the network, rather than closer to 
the sunspot, however still in structures that appear connected 
to the sunspot. Figure [5] shows in the upper panel a frequency 
histogram of the projected velocity determinations. 

The distribution of all measurements peaks around 
20 km s _1 , with a tail extending up to almost 80 km s _1 , re- 
sulting in an average velocity of 27.4 km s _I . There is a 
significant difference between inflow (towards the sunspot) 
and outflow, yielding average velocities of 30.2 km s _1 and 
20.0 km s _1 , respectively. The high- velocity tail of the distri- 
bution is dominated by features moving towards the sunspot, 
while the lower velocities are relatively more abundant for 
features moving in the opposite direction. 

3.1.2. Doppler Velocities 

As the Gaussian fit requires a minimum of three data points, 
this results in fewer successful Doppler velocity determina- 
tions than for the first moment method, reducing the total 
number of measurements by about 18%, compared to the re- 
duction of about 9% with the first moment approach. Never- 
theless, both the distributions of Doppler velocities are very 
similar regardless of the method, their averages falling within 
each-other's measurement error as determined above. On one 
hand, this validates the first moment approach, while showing 
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Fig. 5. — Velocity histograms of the flocculent features with a bin size of 
5 km s _1 . Top panel: Projected velocities obtained from measurement in the 
space-time diagrams. Middle panel: Doppler velocities. Bottom panel: Ab- 
solute velocities obtained by combining the projected and Doppler velocities 
through vector addition. In each panel the solid distribution describes all 
measurements, the dashed one the subsample of features moving towards the 
sunspot, the dotted one those moving away. The vertical dashed lines indi- 
cate the average velocity all considered measurements, while the dotted lines 
in the middle panel show the average Doppler velocity of the positive and 
negative subsamples. 

on the other hand that the influence on the derived statistics of 
the reduced number of valid measurements for the Gaussian 
fit method is small. 

The middle panel of Fig.[5]shows the overall Doppler distri- 
bution, as well as that of the subsamples for features moving 
towards and away from the sunspot. The determined Doppler 
velocities range between -31.4km s _1 and 36.6 km s _1 , with 
an average of -8.2 km s _1 (or -16.5 km s _1 and 14.6km s 
when considering the inflowing and outflowing subsamples, 
separately). The blueshift peak is dominated by features 
moving towards the sunspot, as expected from the orienta- 
tion of the majority of the tracks shown in Fig. Q] with re- 
spect to the slanted line-of-sight. Conversely, features that 



move away from the sunspot along these tracks are red- 
shifted. The average Doppler shift for the blueshifted features 
is -22.6 km s , with the average of the subsamples differing 
by less than about 0.3 km s from the total average. The red- 
shifted features show a double-peaked distribution, centered 
around 10km s _1 and 35km s _1 , respectively, and with an av- 
erage of 21.9km s . As the contribution of the outflowing 
features is mostly to the lower of those two peaks, its average 
is correspondingly lower than that of the inflowing features, 
that contribute almost exclusively to the higher velocity peak. 

3.1.3. Absolute Velocities 

Combining the projected and Doppler component veloci- 
ties produced the distributions in the bottom panel of Fig. [5] 
Both the shapes of the distributions and range in velocities 
are similar to those in the upper panel. The absolute veloci- 
ties found range between 7. 2-8 2. 4 km s _1 , with an average of 
36.5 km s _1 when taking all measurements into account. Con- 
sidering the inflow and outflow subsamples separately, the av- 
erages are 38.8km s _1 and 30.3 km s , respectively. Prop- 
agating the uncertainties determined for the projected and 
Doppler velocities to the absolute velocity, the average error 
in the absolute velocity was found to be 5.9 km s _1 . 

3.1.4. Accelerations 

Figure [6] shows the frequency histogram for the 50 accel- 
eration measurements. Of these, 40 correspond to acceler- 
ations, whereas 10 are of decelerations (of which 8 are of 
features moving towards the sunspot and 2 of features mov- 
ing in the opposite direction). The average acceleration (i.e., 
only including positive values) is 270 m s~ 2 , which lies close 
to the gravitational acceleration in the solar photosphere of 
274 m s~ 2 . However, this coincidence is fortuitous since most 
values are much smaller; the average is boosted by the eight 
outliers with values in excess of 400 m s~ 2 . Note that the mo- 
tion of the flocculent features is likely field-guided, so that 
only the component along the fibrillar path acts as effective 
gravity. 
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Fig. 6. — Histogram of the acceleration measurements. Solid line: all mea- 
surements. Dashed line: features moving towards the sunspot. Dotted line: 
features moving away from the sunspot. The vertical lines indicate the av- 
erage over all measurements {dashed), as well as decelerations and accelera- 
tions (dotted). The bin size is 50m s -2 . 

The average deceleration is much smaller (149 m s~ 2 ), with 
an only slightly higher value of 162 m s~ 2 when considering 
the features moving towards the sunspot and a yet smaller 
value (96 m s~ 2 ) for features moving towards the network. 
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3.1.5. Lifetimes and Occurrence Frequency 

Lifetimes were determined using the start- and endpoints 
of each space-time feature track as in Fig. @] The right-hand 
panel of Fig. Q shows the resulting distribution. The 12.4 s 
lower cutoff results from requiring at least two samplings 
for a velocity measurement. The mean visibility lifetime is 
67.7+8.8 s, with no significant difference between in- and out- 
ward motion. Note that the lifetimes are rather lower limits, 
as the features become obscured by fibrils at line center when 
the Doppler velocity becomes small. 

The occurrence frequency was derived by counting the 
number of measurements at each single time step, resulting 
in the histograms shown in the left-hand panel of Fig. [7] On 
average 68 features can be observed each minute, with an av- 
erage occurrence rate of 47 and 20 features per minute for the 
in- and outflow subsamples, respectively. The occurrence fre- 
quency appears to decrease somewhat during the second half 
of the observations, after a peak of 145 occurrences in one 
minute. 

3.2. Shapes of Flocculent Features 

Figure |8] shows the size distributions of the flocculent fea- 
tures. Out of 236 features with successfully determined 
Doppler velocities, lengths could be determined for 206 fea- 
tures, while 209 widths were obtained. The lower cutoff in 
both distributions is at the telescope diffraction limit. Up- 
per cutoffs were visually set at 20 px (or about 1 Mm) for 
the lengths and at lOpx (equivalent to slightly over 500 km) 
for the widths. On average, the flocculent features measure 
627 km in length and 304 km in width, but lengths range be- 
tween 353-1005 km, while widths of 141-509 km are found. 

We also determined eccentricities for those 186 features 
with both length and width measurements. In general the fea- 
tures are largest in the propagation direction; only 3 features 
have a larger "width" than "length". On average, the eccen- 
tricity of the flocculent features is 0.83+0.05, with the same 
value for the average for the larger prolate sample of 183 fea- 
tures, while the average for the oblate sample of 3 features is 
much smaller (0.27+0.05). 

The morphology of the flocculent features changes only 
slightly as they progress along their trajectories. Some fea- 
tures appear to get stretched in length as they move towards 
the sunspot, while others appear to be compressed. Slightly 
more features are found to expand rather than shrink in length 
as they propagate. 

3.3. Doppler Shift of the Line Core 

Figure|9]shows the time-averaged result of the core Doppler 
shift measurements. The average values lie between about 
-5.7 km s _1 and 13.5 km s _1 , where the minus sign (and blue 
color) indicates motion towards the observer. The average 
Doppler velocity over the full field-of-view is 2.7 km s _1 , di- 
rected away from the observer. 

There is a conspicuous asymmetry near the sunspot be- 
tween blueshifted patches on the limb-ward side and red- 
shifted ones on the disk-center side. These illustrate gen- 
eral inflow on the order of a few kilometers per second, with 
slightly higher velocities of up to 10 km s in the redshifted 
patches and up to -5 km s _1 in the blueshifted ones. The 
blueshifted limb-ward patches generally coincide well with 
the ends of the trajectories of the flocculent flows, as well 
as with the superpenumbral fibrils observed in Ha line cen- 
ter (cf. Fig. [TJ. On the disk-center side such correlation is 



less strong. The temporal evolution of line core Doppler shift 
maps shows the blueshifted patches on the limb-ward side of 
the sunspot to be relatively stable, while the redshifted patches 
on the disk-center side are more variable, being replaced by 
blueshifts for several minutes at a time before returning to a 
redshift. 

Space-time diagrams extracted from the Doppler map cube 
along the same flocculent flow trajectories show only occa- 
sionally enhancements that are co-spatial with the dark streaks 
observed in the space-time diagrams from the regular inten- 
sity data. In those cases for which corresponding enhance- 
ments are found, the related (absolute) line core Doppler shifts 
are typically on the order of 5-15 km s _1 . Also, the sign of the 
core Doppler shift correlates well with that of the flocculent 
feature Doppler velocity. 

4. DISCUSSION 

Our key finding is the presence of high-speed flows towards 
and away from the sunspot that are of intermittent, flocculent 
nature rather than continuous. They are observed in the He 
wings and appear to follow the fibrillar chromospheric struc- 
tures seen in the Ho- core. Summarizing their properties: 

1. We observe dark flocculent features with an average 
length of 627+44 km and width of 304+30km propa- 
gating through the chromospheric canopy. 

2. An average absolute velocity of 36.5+5.9 km s is 
found, but velocities in the range of 5-80 km s are 
also measured. Features moving away from the sunspot 
propagate at a significantly lower average velocity of 
30.3 km s , against an average inflow at 38.8 km s _1 . 

3. On average, flocculent features can be observed for 
67.7+8.8 s and typically 68 features are visible each 
minute. The actual lifetimes may be considerably 
larger, as the features typically (dis)appear in mid-flight 
or before reaching the photosphere as observed in the 
far wings. 

4. Though practically co-spatial, the flocculent flows ap- 
pear to be a phenomenon distinct from the "classical", 
lower-velocity inverse Evershed effect. 

In this section we first discus s the se flows with resp ect to 
the inverse Evershed effect (Sect. 14. il l, coronal rain (Sec t.l4.2T >. 
and then discuss the possible driving mechanism (Sect. l4~3l >. 

4. 1 . Flocculent Flows and Inverse Evershed Effect 

The general picture of inflow in the chromospheric super- 
penumbra is confirmed by the results obtained in this study. 
Not only is it observed as an inflow at line -of-s ight veloci- 
ties of a few kilometers per second (cf. Sect. 13.3b , but also in 
the form of flocculent features moving at considerably higher 
velocities. Although we find slightly higher average values 
of up to 15 km s _1 , the line core Doppler shift velocities are 
similar to the inverse Ev ershed effect ve locities obtained in 
earlier studies. Already ISt. John! (11913b found inverse Ev- 
ershed velocities on the order of about 3 km s" 1 on aver- 
age an d similar velocities were obtained by lBones & Maltbyl 
(119781) . who measured an average radia l line-of-sight ve- 
locity of 1.8 km s" 1 from Mgbj spectra. iDere et all (119901) 
studied Hq±0.3 A velocity maps and foun d velocities of 
2.6km s _I , while [Alissandra kis et al.l {1988) obtained chro- 
mospheric line-of-sight velocities (averaged over the posi- 
tion angles considered) between +2 km s (i.e., upflow) and 
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Fig. 7. — Occurrence frequency (left panel, 1-min bins) and lifetimes (right panel, 12.4 s bins). Solid line: all measurements. Dashed line: features moving 
towards the sunspot. Dotted line: features moving away from the sunspot. The vertical dashed line in the right-hand panel indicates the average lifetime. 
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i, — Size histograms showing flocculent feature lengths (left panel, bin size of 100 km) and widths (right panel, bin size of 50 km). 

-4 km s _1 , and lHaugenl (fT969s!) found an average velocity vec- 
tor peaking at the penumbral rim at 6.8 km s _1 . The veloc- 
ities measured for the flocculent features, a re more on the 
10 ord er of the chromospheric inflow velocities Maltby] (1 1975b 

and iBeckersI (1 1962b found in their filtergram studies (veloc- 
ities of 5-15kms _I and 40-50 km s , respectively), but 
^ are also comparable to the velocities of some of the single 
5 "I absorbing elem ents measured in the spectroscopic study by 
"=< iHaugenl ( I 19691) . who found velocities ranging between -35 
f and +50 km s _1 (where positive values indicated downflows). 
s? We suggest nevertheless that most of the older studies, ham- 

pered by much lower spatial resolution than reached here and 
considering their morphological description of the inflows, 

concern the more gentle flow patterns that we observe as 
patches of enhanced line core Doppler shift seen in Fig. [9] The 
flocculent features described here are in fact not observable at 
-5 line center where they seem to be obscured by overlying fib- 

rils (cf. Fig. |2). The flow channels of the inverse Evershed 
effect and the flocculent flows appear to be distinct and are 
thus not as such the same flows. 
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Fig. 9. — Time-averaged core Doppler shift map of the full field-of-view. 
Red: positive Doppler velocity (redshift). Blue: negative Doppler velocity 
(blueshift). Dotted: the same tracks as in Fig.[T] Solid: contours for the 
umbra-penumbra boundary and the outer penumbra boundary at Ha -1.1 A. 



4.2. Flocculent Flows and Coronal Rain 

Morphologically, the flocculent flows are rem- 
iniscent of coronal rain. Recent work by 

lAntolin & Rouppe van der Voorll ( 120121) and in particu- 
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lar lAntolin et al.l (120121) shows this is not only the case 
qualitatively, but also quantitatively. Both studies consider 
similar Ha imaging spectroscopy data obtained with CRISP 
as for this study. lAntolin & Rouppe van der Voortl (120 121) 
suggest that the coronal rain phenomenon might be a 
common one, with single condensations attaining average 
absolute velocities on the order of 70km s and typical 
sizes of 0.74 Mm in length and 0.31 Mm in width. Similar 
dy namics and morpho logical quantities have been found 
by lAntolin et al.l (120121) for on-disk cases of coronal rain, 
including also in the currently considered data set. As these 
similarities may lead to difficulties in distinguishing both 
phenomena, as well as being indicative of a possibly related 
driving mechanism, we discuss this phenomenon in more 
depth. 

Typical velocities for the coronal phenomenon are on the 
order of a few tens to a fe w hundred kilom eters per sec- 
ond, as for instance found bvlSchriiverl (1200 1[) ( studying qui- 
escent coronal loops) and lO'Shea et al.l (120071) . who found 
Doppler shifts of 50- 100 km s _1 and velocities of at least 
up to 100 km s _1 from He i and O iv data obtai ned with the 
Coron al Di agnostics Spect r ometer on SOHO. Muller et al. 
( 120051) and iDe Groof et al.l (120051) obtained velocities in a 
similar range (30 - 120km s~') from observations, although 
the simulations of Muller et al.J ( 120051) yielded typical veloc- 
ities of about 50 km s -1 , while also indicating that subse- 
quently downflowing blobs may travel at increasingly higher 
velocities of up to 230 km s due to evacuation of the loop 
by preceding condensations. 

The flocculent features ob served in this study are somewhat 
smalle r tha n those found by lAntolin & Rouppe van der Voortl 
(120121) and lAntolin et al.l (12012b . however mostly so in length, 
and move at much lower average velocities than coronal rain. 
The lower velocities may result from the flocculent flow chan- 
nels being much shorter than coronal loops, as the traced tra- 
jectories range between 3-14 Mm. iMiiller et al.l (120031) re- 
ported on simulations of 10 Mm coronal loops, but found 
propagation velocities on the order of only a few kilometers 
per second for the condensations that formed. 

A further, important, difference is that flocculent features 
appear to constitute a footpoint-to-footpoint, rather than apex- 
to-footpoint, flow. Even though no flocculent feature could be 
traced all the way from one footpoint to the other in this data 
set, the (re)appearance of the features close to the fibril apices 
moving at considerable velocities would argue against them 
forming through catastrophic cooling alone. In that respect, 
the highe r flocculent flow vel ocities compared to the veloci- 
ties from IMiiller et al.l d2003l) can be understood, as coronal 
rain accelerates from zero velocity as it starts precipitating. 

Yet, even though the exact mechanism causing coronal rain 
may not be applicable to the flocculent flow phenomenon, 
a similar morphology may be nonetheless indicative of a 
thermal non-equilibrium process leading to the observed 
condensations. Comparison wi t h the results obtained by 
lAntolin & Rouppe van der Voortl (120121) shows yet another 
striking similarity between both phenomena: in a majority of 
the flocculent flow cases adjacent fibrils display blobs prop- 
agating in the same direction at the same time, which may 
point at co-evolution of the conditions in neighboring fibrils. 

4.3. Driving Mechanism 

Distinguishing wave phenomena from actual mass flows 
is a recurring challenge, as what is observed as floccu- 
lent features may actually be gas compressed by prop- 



agating wave fronts, their velocities reflecting the phase 
speed of the waves, rather than actual mass motion. In- 
deed, similar features that have been ob served before by 
IVan Noort & Rouppe van der Voortl (12006ft . who studied Ha 
line center data (i.e., no further spectral information), were 
explained as MHD waves. The qualitative morphology and 
behavior of the "moving features in loops" is practically iden- 
tical to the flocculent flows described here and the differences 
with the projected velocities they report (40-120 km s~') may 
well be attributed to projection effects. With typical sizes of 
up to 275 km, the features they analyzed are notably smaller, 
however. Also, they report on bright blobs at line center, 
whereas the flocculent flows can only be seen in the wings 
of Ha, despite the viewing angles being comparable for both 
data sets. This may, however, be due to an inherent differ- 
ence in fibril topology or by a different orientatio n of those 
fibrils to the line-of-sight. A more recent study by iLin et"aT] 
(1201 ll) . covering in part the same data sets, expands on the 
earlier results by IVan Noort & Rouppe van der Voortl (120061) 
by including CRISP data. They show that the bright blob-like 
features typically have small Doppler shifts (+5-10km s~'), 
last for up to about 150 s, and slightly increase in length as 
they propagate. 

Features with similar morphology and behavior have also 
been found in other st udies, albeit propagating at lower ve- 
locities. For instance, Sanc hez- Andrade Nuno et ail d2008l) 
found chromospheric magneto-acoustic waves in Ha ob- 
servations. These waves manifest themselves mainly as 
(trains of) blobs with average dimensions as small as 725 km 
(1 arcsec) in length and 360 km (0.5arcsec) in width, prop- 
agating at phase velocities between 12^-2kms _1 . Both 
the sizes and the phenomenology are similar to what we 
observe for the flocculent flows. However, in contrast to 
th ose (but in common with some of t he fea tures described 
by I Van Noort & Rouppe van der Voorq (|2006l) X the blobs ob- 
served by ISanchez-Andrade Nuno et all show displacement 
perpendicular to their propagation direction. In a more re- 
cent study of Ca n 8542 A data obtained with SST/CR1SP 
IWedemeyer-Bohml (1201 Oh reported on (among other things) 
bright blobs propagating at projected velocities of 4-8 km s _1 
along neighboring fibrils that he explained as a wave manifes- 
tation. 

Considering this scenario further, we note that typ- 
ical chromospheric sound speeds are on the order of 
1 0-15 km s 1 (e.g.. lAlissandrakis et al.l 1988b ITsiropoulal 
12000b ISanchez-Andrade Nuno et all 120081: ILin et alJ 1201 W . 
while Alfven speeds may range between 175-600 km s _I 
for magnetic field strengths b etween 30-100G (cf. 
ISanchez-Andrade Nuno et all 120081 and references therein) 
and a reasonable chro mospheric density of 2xl0~ 13 g cm 
(cf. ITsiropoulal 120001 and references therein). Although 
subsonic flocculent features have been found, most exceed 
the chromospheric sound speed by far and hence the floc- 
culent flows cannot be attributed exclusively to acoustic 
waves. A magneto-acoustic wave scenario may allow for 
much higher velocities (even compa rable to flocculen t flow 
velocities), as for instance shown by iGiovanellil (119751) . who 
measured phase velocities on the order of 70 km s in Ha 
observations of quiet chromosphere and a sunspot. Assuming 
typical chromospheric number densities and a magnetic 
field strength of 10G, these velocities could be reconciled 
with an Alfven wave scenario. However, one thus needs to 
assume rather low field strengths (i.e., below 6 G) in order 
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to obtain an Alfven velocity as low as the average flocculent 
flow velocity of 36 km s -1 . In addition, the high Doppler 
velocities of the flocculent features in combination with the 
changes in the Doppler component velocities as the features 
propagate along their curved paths, would seem to argue for 
actual mass motion, rather than waves. 

An alternative driving mechanism could be found in the 
siphon flow model. Siphon flows have been sugge sted be- 
fore in the context of the inverse Evershed effect ( cf. iHaugenl 
19691: iMaltbvl [1971 lAlissandrakis et al.l [1981 iDere et alj 
199d Boerner & Kned Il992h iKjeldseth-Moe et all 11993b 
Montesinos & Thomasll 19971) . In this model, a flow is driven 
by a gas pressure inequality between the footpoints of a mag- 
netic flux tube, which in turn is caused by a difference in mag- 
netic field concentration (and hence strength) at the respective 
footpoints. A higher magnetic field strength at one footpoint 
will result in a lower gas pressure at that same footpoint and 
vice versa. As the magnetic field strength of a sunspot will 
be larger than in the surrounding plage, a flow directed to- 
wards the sunspot would naturally occur in this model. An- 
alytical solutions predict both sub- and supersonic flows, de- 
pending on the parameters of the loop structure, however su- 
personic flows in the critical solution (i.e., where the flow 
has been accelerated from sub- to supersonic) will always re- 
turn to a subsonic state through the fo rmation of a shock (cf . 
ICargill& Priestl 1 19801: iThomasI [19881 Since the flocculent 
flows appear to be footpoint-to-footpoint flows, the siphon 
flow model may offer a viable way to explain them. For most 
cases a supersonic flow would be required, at the very least for 
that part of the fibrils where the flocculent features are visible. 

A number of problems arise, however. The siphon flow 
model may explain the general characteristic of an inward 
flow towards the sunspot, but it does not explain why floc- 
culent features could be formed. Although shocks will occur 
in the transition from supersonic to subsonic flows and have 
been observed to a larg e r or smaller degre e of c onfidence (cf. 
Degenhardt et al.1119931: lUitenbroek et ai1l2006t iBethge et al.l 
20121) . there is no evidence that these shocks cause propa- 
gating condensations and would rather show up as brighten- 
ings due to material compression and the consequent rise in 
temperature. A few cases of decelerating flocculent features 
were found in this study, but the return to subcritical veloc- 
ities th rough a shock was not evident from the data. How- 
ever, as lBethge et al.l (120121) already pointed out for their case, 
these shocks could be occurring at lower altitudes that are not 
sampled by the Ha line. Alternatively, the shocks could be 
obscured from view by overlying fibrils, in a similar way to 
flocculent features closer to Ha line center. Also, we observe 
flows towards the plage, i.e., in the direction opposite to what 
the siphon flow model predicts. Although in most cases these 
features are seen close to the plage (footpoints) only and could 
thus be propagating in fibrils that are actually not connected 
to the sunspot (allowing possibly low enough field strength 
for oppositely directed flows), this d oes not hold true for al l 
cases. Nonetheless, as sugg ested by |Cargill& Priestl dl980l) . 
an oppositely directed flow may occur in a fibril with an es- 
tablished siphon flow as a result of a change in pressure differ- 
ence. Evidence for the exi stence of such a counterflow has fo r 
instance been provided by iGuglielmino & ZuccarelTol (1201 ll) . 
in a recent study of a pore and its surroundings. 

Although consid ering much larger scal es, recent simula- 
tions described by IZacharias et alj d201 ll) may also offer a 



possible driving mechanism. The simulations they present 
show the emergence of a blob-like feature from the transition 
region, resulting from the build-up of pressure (in turn caused 
by a heating event at one footpoint) propelling the feature 
further through the corona along a magnetic loop structure 
before falling down again. The feature described has much 
larger dimensions than typical flocculent features or coronal 
rain and only occurred once during a one hour simulation 
run. A similar heating mechanism might nevertheless work 
on smaller scales too, producing blobs of smaller sizes prop- 
agating through the chromosphere and/or transition region at 
possibly shorter time intervals. 

5. CONCLUDING REMARKS 

In this paper we have analyzed flocculent flows occurring 
along the chromospheric canopy as observed at high spatial, 
temporal and spectral resolution in Ha data obtained with the 
CRISP instrument at the Swedish 1-m Solar Telescope. Al- 
though we have been able to obtain statistical properties of 
a considerable sample of flocculent features, we can at this 
point merely speculate about the mechanism driving the floc- 
culent flows. The results obtained, both qualitatively and 
quantitatively, seem to argue against a purely wave-related 
phenomenon. Rather, the flocculent features are more likely 
to occur as part of an already present (supersonic) siphon flow, 
which could develop its flocculence as a result of, for instance, 
recurring heating events leading to the formation of conden- 
sa tions (in analogy to the larger scale phenomenon described 
by IZacharias et alJ(l201 lb ), or waves interacting with the flow- 
carrying fibrils, causing the condensations to form as the fibril 
is pushed out of thermal equilibrium. In either case, these con- 
densations would subsequently be carried by the siphon flow 
at supersonic velocities (where they are observed) before be- 
ing decelerating through a shock to subsonic velocities upon 
reaching the photosphere. 

Given that this phenomenon appears to be ubiquitous and 
not solely related to the superpenumbra, a study of multiple 
data sets (including also more quiet regions) and with different 
diagnostics would not be out of place. In addition, where ob- 
servations may constrain the properties of the flocculent fea- 
tures as well as the necessary environmental conditions for 
them to occur, simulations would be very helpful in further 
uncovering their driving mechanism. 
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ANALYSIS TOOLS: INTRODUCING CRISPEX AND TANAT 

Sophisticated instruments such as the CRISP instrument at the SST deliver large data sets that are often formatted as multi- 
dimensional data cubes. These are nontrivial to browse and analyze, and require therefore a tool with which one may quickly 
explore the data obtained for interesting features and which, preferably at the same time, allows for more in-depth analysis of the 
data set of interest. A widget based tool, programmed in the Interactive Data Language (IDL), was developed with that goal in 
mind and driven by the particular challenges faced as a result of the extensive observational possibilities offered by CRISP. The 
CRisp SPectral EXplorer (CRISPEXfl can currently handle: 



• single spectral scans; 



• simple 3D temporal cubes (i.e., a time series of intensity images); 



• spectrotemporal cubes (i.e., cubes with both spectral and temporal information); and 

• Stokes data cubes (i.e., spectrotemporal cubes including also multiple Stokes parameters). 

As the typical size of CRISP data cubes exceeds the available random access memory on regular laptop or desktop computers, 
the data are read in to an associated variable and hence CRISPEX need neither read nor keep a full data cube in memory, thereby 
enabling quick browsing in the spatial, temporal and spectral domain. This swiftness comes at a cost, however. In order to 
have full spectral browsing functionality, a spectral data cube (which is nothing more than a reordered image data cube) must be 
supplied as well, thus requiring double the amount of disk space for one data set. 

Though originally designed and intended to handle CRISP data, CRISPEX can cope with any (observational or synthetic) 
data, provided it has been formatted in a certain way. Indeed, the tool has already been successfully used in several studies 
based on CRISP data dRouppe van der Voort eFa l. 2009; Wed emever-Bohm & Rouppe van der Voorlll2009t lAntolin et alJl20Tot 
IWatanabe et"aT]l20lU iLin et alJl201 ll) . as well asHinode/SOT data (lAntolin & Verwichtdl201 ll) . 

Visualization and Analysis Options 
Depending on the input data, different visualization and browsing options are available to the user. One can simultaneously: 



• play forward and backward in time; 



• step through a spectral profile; 



• switch between Stokes parameters; 



• blink between image frames, stepping both in temporal and spectral dimension. 

Data feedback is given to the user through multiple screens, of which Fig.[lO]shows an example. Following CRISP's capability to 
run multiple line programs (i.e., observing programs where time series in more than one spectral line are obtained), CRISPEX can 
display two data sets at a time. Such reference data (e.g., a reference magnetogram) need not have the same spectral dimension, 
but are currently only constrained to having the same spatial dimensions and either consist of a single frame or the same amount 
of temporal frames as the main data cube. Additionally, CRISPEX allows for simple image and spectrum saving (with the option 
to include or exclude overlays, such as the cursor or drawn paths), which can be used for simple movie production. 

CRISPEX also offers a number of analysis options, but one of its first and most developed features is the extraction of 
space-time diagrams along user-defined (curved) paths. These paths can (as they are in this study), but need not necessar- 
ily be, manually defined. One m ay also supply a file with detections to extract space-time diagrams along, as was done by 
iRouppe van der Voort et al.l ((2009) in their spectral analysis of on-disk counterparts of type II spicules. Space-time diagrams can 
be extracted and saved for a single wavelength position or for a given spectral range (which may be the entire spectral scan), 
as well as for a restricted temporal range. Another analysis option is given by the spectral slit functionality, which extracts the 
spectrum for each pixel along a line centered on the cursor. The slit can be adjusted both in length and in rotation angle, thereby 
allowing the coverage of practically any feature of interest. In particular for those cases where a profile scan is available, the 
in-program Doppler image functionality may be of use to browse the Doppler image constructed from the subtraction of the red 
from the blue wing. Though this feature may slow down the overall playback, depending on the memory available on the machine 
used, one escapes the need to produce a separate Doppler cube. Lastly, the functionality to display intensity-time diagrams may 
be particularly useful when cubes do not consist of profile scans but rather of different spectral diagnostics. This functional- 
ity could be used to show, for example, subsequent detection of a phenomenon in different diagnostics of multiple SDO/AIA 
channels. 



Both the actual distribution and more detailed descriptions of the func- | http7/bit.lv/crispex I 
tionality and options of CRISPEX and its auxiliary routines can be found at 
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CRISPEX: Parameters overview 



trar x: 368 Wsvelengtri index: 18 Stokes: I ^ 

Cursor y: 637 Wauelengtk: Tlrae index: 27 

roomfactor: 1 Dcprler oolocitu (kr/s): 3.1 .33 netusi time (si: iS7.4u 



Fig. 10. — Default layout of windows loaded by CRISPEX, showing clockwise from the top left: (1) the intensity image, (2) the spectrum-time diagram, (3) the 
local spectrum, (4) the control panel and (5) the parameters overview window. 

Auxiliary Programs 

Although the CRISPEX distribution comes with several auxiliary routines, one program in particular deserves some further 
introduction, since it offers important analysis options complementary to those in CRISPEX. The Timeslice ANAlysis Tool 
(TANAT) facilitates the analysis of space-time diagrams extracted with CRISPEX, as it enables the determination of projected 
velocities and accelerations of features from these diagrams. Such measurements can subsequently be saved to a file, with the 
option to add a (custom) flag for each measurement as it is saved, and adding by default a unique identifying number which 
enables the correct overlay of saved measurements on the corresponding space-time diagram. In addition, different visualization 
options allow switching between or combination of space-time diagrams for different positions within the spectral line observed, 
thereby enabling - for instance - the comparison of a feature's signal (or lack thereof) in the red and blue wing of a spectral line. 
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